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1. Introduction 

In this project, we perform first steps towards extracting the properties of the £?o(980) and k 
resonances from lattice QCD. 

The mass and width of an infinite volume resonance can theoretically be extracted from Eu- 
clidean lattice QCD by studying the corresponding energy levels as a function of the lattice size L 
with an extension of Liischer's finite volume method [fil 51 0, ^]. 

However, such an analysis is challenging because several energy levels have to be extracted 
from the lattice simulation with high accuracy. Therefore, we follow an exploratory approach: we 
perform a study of correlator matrices of four-quark operators in order to investigate whether bound 
four-quark states are observed in the ao(980) or K channels. Recently, in such a set-up in Nf = 2 
QCD, hints at bound tetraquark states in the a and K channels have been identified [g]. 

We compute the correlator matrices in a mixed action [J7|, |8|] analysis of four-quark states on 
the lattice ignoring disconnected contributions. The goal is to answer the question whether there 
is a bound tetraquark or molecule state in the ao(9SO) or K channel. Identifying such a state in 
addition to the free scattering states could be a hint on the nature of the corresponding resonance. 
Additionally, the overlap to different interpolating operators could be studied. 

The lattice study is based on Nf ■ = 2 + 1 + 1 gauge configurations generated by the European 
Twisted Mass (ETM) Collaboration. Details on the generation of the ensembles can be found in 
Ref. [^]. The analysis is performed on four ensembles with a lattice spacing of a « 0.086 fm and 
pion masses between m n + 280MeV and m n + 450MeV, see Table [|. 

For the computation of observables we use a twisted mass discretization for the valence s 
quarks, which is different from the sea s quarks to avoid the problem of mixing between s and c 
quarks, for details cf. Refs. 0>§]. 

2. Interpolating operators and extraction of energy levels 

To perform a four-quark analysis of ao(980), we compute correlators C,- 7 (?) =< 0,(?)C?j(O) > 
with interpolating operators 

°* = E £ {d{x)^s(x)){s{x)Tfu{x)) fori = 1,2,3, (2.1) 
Oi = £ [d(x)Yif{x)] a [s T {x)T i u{x)] a for i = 4,5, (2.2) 

X 

6 = £ (s(x)y 5 s(x)) (d(x)y$u(x)) , (2.3) 

x 

with the y-matrices F\ = y 5 , F% = y M , = 7^75, T4 = C75, F 5 =C, C = 7072. Due to the omis- 
sion of disconnected diagrams the artificial pseudoscalar sy^s state, which we call T] s in the follow- 
ing, will be relevant for the calculation. To ensure the correct identification of the Tj s 7r state, we 
explicitly included the operator Of,. 
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ensemble 


P 




a\i a 


a^s 


L/a 


Aconf 


smearing 


A30.32 


1.90 


0.0030 


0.150 


0.190 


32 


672 


APE 


A40.20 


1.90 


0.0040 


0.150 


0.190 


20 


500 


none 


A40.24 


1.90 


0.0040 


0.150 


0.190 


24 


1259 


APE 


A80.24 


1.90 


0.0080 


0.150 


0.190 


24 


1225 


APE 



Table 1: The input parameters of the ensembles used in this project, the number of configurations and the 
smearing type employed. 

For further clarification of the nature of the examined states, the two-meson operators 

On = £ {d(x)y 5 s(x)) £ (s(y)Y 5 u(y)) , (2.4) 

x y 

8 = £ (s(x)y 5 s(x)) £ {d(y)Ysu(y)) (2.5) 

x y 

are studied as well (on a single ensemble). 

In order to extract energy levels from the correlator matrices, we solve the generalized eigen- 



value problem (GEVP) [Q, [TO], TJJJ . However, due to the presence of pairs of pseudoscalar me- 
son states \M\ > and \M-2 > coupling to our operators according to < M\\Oj\M2 >^ problems 
arise: one of the two mesons can travel forward the other backward in time, severely complicat- 
ing the analysis of our correlators. To avoid these problems, we restrict the analysis to t < T/4 
and T — t < T/4, because the single-meson contributions are not relevant for sufficiently small t 
(or T — t). For a more detailed discussion we refer to an upcoming publication and to the Refs. 



12, 13]. Excited state contributions are taken into account by performing two-mass fits to the 



eigenvalues of interest. 
3. Results 

3.1 ao(980): four-quark and two-particle operators, a single ensemble 

We start by discussing ao(980) (I(J PC ) = 1(0 ++ )) results obtained using ensemble A40.20 
(cf. Table [I]). This ensemble with rather small spatial extent (L 1.72fm) is particularly suited to 
distinguish two-particle states with relative momentum from states with two particles at rest and 
from possibly existing ao(980) four-quark states (the former have a rather large energy because one 
quantum of momentum is p m i n = 2n/L 720MeV). 

Figure [j]a shows effective mass plots from a 2 x 2 correlation matrix with a KK molecule 



operator 0\ (see equation (|2JJ)) and a diquark-antidiquark operator O4 ( |2.2[ ). The corresponding 
two plateaus are around 1000 MeV and, therefore, consistent both with a possibly existing ao(980) 
four-quark state and with two-particle K + K and r) s + 7i states, where both particles are at rest. 

Increasing this correlation matrix to 4 x 4 by adding two-particle K + K and T] s + % opera- 
tors (equations (2.4) and ( |2.5[ )) yields the effective mass plot shown in Figure |l|b. The same two 



low-lying states are resolved, however, with significantly better quality. Two additional states are 
observed, whose plateaus are around 1500 MeV — 2000 MeV. From this 4x4 analysis we conclude 
the following: 
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Figure 1: ao(980) sector, A40.20 ensemble, (a) Effective masses as functions of the temporal separation, 
2x2 correlation matrix (local operators: KK molecule, diquark-antidiquark). Horizontal lines indicate 
the expected two-particle K + K and T} s + % energy levels, (b) 4 x 4 correlation matrix (local operators: 
KK molecule, diquark-antidiquark, two-particle K + K, two-particle rj + it), (c), (d) Squared eigenvector 
components of the two low-lying states from (b) as functions of the temporal separation. 



1. We do not observe a third low-lying state around 1000 MeV, even though we provide op- 
erators, which are of four-quark type as well as of two-particle type. This suggests that the two 
low-lying states are the expected two-particle K + K and f] s + 7l states, while no additional stable 
^0(980) four-quark state does exist in the A40.20 ensemble. 

2. The effective masses of the two low-lying states are of much better quality in Figure |I|b 
than in Figure [l]a. We attribute this to the two-particle K + K and r\ s + % operators, which presum- 
ably create larger overlap to those states than the four-quark operators. This in turn confirms the 
interpretation of the two low-lying states as two-particle states. 

3. To investigate the overlap in a more quantitative way, we show the squared eigenvector 
components of the two low-lying states in Figure [T]c and Figure |l|d (cf. Ref. [||] for a more detailed 
discussion of such eigenvector components). Clearly, the lowest state is of r/ s + % type, whereas the 
second lowest state is of K+K type. On the other hand, the two four-quark operators are essentially 
irrelevant for resolving those states. These eigenvector plots give additional strong support of the 
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above interpretation of the two low lying states as two-particle states. 

4. The estimated energy levels of two-particle excitations with one relative quantum of momen- 
tum are consistent with the effective mass plateaus of the second and third excitation in Figure |T]b. 



Figure [T]a and Figure [T]b also demonstrate that two-particle states can be resolved by four-quark 
operators, i.e. two-particle operators are not necessarily needed, to extract the full spectrum. Since 
we are mainly interested in possibly existing states with a strong four-quark component, we restrict 
the correlation matrices computed for other ensembles to four-quark operators. 

3.2 «o (980) : four-quark operators, many ensembles 

We have analyzed the three additional ensembles listed in Table [I] with respect to (20(980) in 
a similar way as explained in the previous subsection. The main difference is that this time we 
exclusively use four-quark operators, but no two-particle operators. To be able to resolve more 
than two low-lying states, we use the operators 0\ to Of,. 

An effective mass plot for the A30.32 ensemble (cf. Table [j]) is shown in Figure ^| together 
with the expected energy levels of the relevant two-particle states and the masses extracted by fits. 
The effective mass plots for the other ensembles will be shown in an upcoming publication. 

On a qualitative level our findings agree for all ensembles, i.e. are as reported in the previous 
subsection: there are always two low-lying states, whose masses are consistent with the expected 
masses of the two-particle K+K and r\ s + % states; higher excitations (the third, forth, etc. extracted 
state) are in all cases significantly heavier and consistent with two-particle excitations with one 
relative quantum of momentum. 
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Figure 2: Effective masses as functions of the temporal separation for the ao(980) and K sector for the 
A30.32 ensemble. Horizontal lines indicate the expected two-particle K + K and rj s + % energy levels. 
Additionally, the masses extracted from fits are included. 



3.3 fc: four-quark operators, many ensembles 

The analysis for the K sector (I(J P ) = l/2(0 + )) closely parallels the analysis of the ao(980) 
sector presented above. We consider 5x5 correlation matrices containing a K% molecule operator 
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analogue to 0\ and further operators corresponding to O2 to O5 with the appropriate valence quark 
content, so that our operators are essentially identical to those considered in Ref. 

In twisted mass lattice QCD the isospin / is not a quantum number. Therefore, it is not suffi- 
cient to only resolve 1=1/2 two-particle K + % states. One has to take into account also mixing 
with 1 = 3/2 two-particle K + % states, i.e. it is necessary to resolve these low-lying two-particle 
states at the same time (details will be discussed in an upcoming publication). 

An effective mass plot for the A30.32 ensemble is shown in Figure || together with the expected 
energy levels of two-particle K + % states and the masses extracted by fits. While effective mass 
plateaus are consistent with these expected two-particle energy levels, there is no indication of 
any additional low lying state, i.e. of a possibly existing bound four-quark fc state. While this is 
suggested by experimental data, it contradicts the findings of the similar recent lattice study of fc 



4. Summary and Outlook 

We computed the low-lying spectrum in the ao(980) and k sectors by employing trial states 
designed to have a substantial overlap with both two-particle and possibly existing tetraquark states. 
With our ensembles, we did not see additional states beside those that can be identified with the 
expected two-particle spectrum. The next states appear roughly consistent with excitations of the 
first quantum of momentum (2n/L) on top of those thresholds. This is somewhat difficult to rec- 
oncile with the additional state found in Ref. ^ in the K channel, despite the rather similar lattice 
setups. 

We find that the low lying spectrum has essentially exclusively overlap to two-particle trial 
states. This suggests that the states that we see are, indeed, the expected two-particles states at the 
threshold and not tightly bound states either of molecular type or diquark-antidiquark type. 

These conclusions can be strengthened by studying more volumes, by introducing twisted 



boundary conditions [ ]14[ ] and by studying further trial states of different type. As for the latter, 
it will be crucial to combine four quarks with traditional quark-antiquark operators including dis- 
connected diagrams. As for the volume dependence, we plan to use the finite volume formulae of 
Liischer [|], |||, |5[ f|, ||] and their extensions to multiple channels developed in Refs. [15, 16, 1^, |l8| ]. 
At present, our limited number of volumes is insufficient for such an analysis. Corresponding 
computations are in progress. 
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